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ABSTRACT 

X-  \ 

X  Problems  connected  with  the  estimation  of  the  ratio  of  the  means  of  a 
finite  bivariate  population  have  been  considered  in  this  report.  The  usual  es¬ 
timator  of  the  ratio,  based  on  the  means  of  a  bivariate  simple  random  sample 
drawn  without  replacement  (s.r.s.  (w.o.r.)),  has  been  compared  with  estimators 
based  on  alternative  double  sampling  design.  Under  this  design  a  very  large 
s.r.s.  (w.o.r.)  is  drawn  for  measuring  only  one  of  the  variables  and  a  subsample 
(s.r.r.  (w.o.r.))  is  drawn  out  of  the  first  phase  units  for  measuring  the  other 
variable.  Efficiency  and  bias  comparisons  have  been  made  by  subjecting  each  of 
the  competitors  to  the  same  budgetary  constraint.  It  turns  out  that  deviation 

f'  r' 

from  the  usual  set-up  sometimes  leads  to  better  sampling  'strategies  . 
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1.  INTRODUCTION 


An  inference  problem,  not  unconmon  in  practice,  is  the  estimation  of  the 
ratio  R-Y/X  of  two  means  of  a  finite  bivariate  population;  [(X^Y^):  i*  1,...,N] 
where  and  are  measurements  of  two  characteristics  X  and  Y  for  a  member  of 
the  population  identified  by  the  integer  i.  Addressing  to  this  problem,  we  have 

A  ^  A 

confined  ourselves  to  estimators  of  the  form  R  ■  Y/X  based  on  some  suitably  chosen 
random  sampling  scheme.  A  commonly  used  Strategy '  (see  Des  Raj.il963J)  is  to 
draw  a  bivariate  simple  random  sample  without  replacement  (s.r.s.  (w.o.r.)), 

A  A 

(x,  ,y, (x  ,y  )  and  to  choose  X«x  and  Y-y,  the  sample  means.  We  propose 
linn  • 

an  alternative  strategy  in  which  the  units  are  drawn  by  'Double  Sampling'  scheme 
(see  Cochran  11977],  pp.  327-355,  Sukhatme  Il954l,  Murthy  [1977]). 

We  have  made  comparisons  by  subjecting  the  competitors  to  a  common  constraint. 
Ensuing  conclusions  are  relevant  to  situations  in  which  the  availability  of  re¬ 
sources  Ce.g.  budget  or  time)  is  a  dominant  factor,  it  is  a  given  constant  B. 

Let  qQ  be  the  overhead  per-unit  cost  of  selecting  a  unit  from  a  given  popula¬ 
tion  and,  q  and  q  be  the  per-unit  costs  of  measuring  respectively  the  charac- 
x  y 

teristics  X  and  Y  for  the  selected  unit.  The  sample  size  for  the  usual  estimator 

A 

R*y/x  is  then  determined  by 

n<q0+qx  +  qy)  “B*  (1,1) 

In  many  practical  situations  q^  and  q^  differ  significantly,  sometimes  to  an 
extent  that  a  sacrifice  of  one  observation  on  the  costlier  component  may  allow 
considerably  large  numher  of  observations  on  the  cheaper  component.  In  situa¬ 
tions  of  this  kind  it  is  possible  to  obtain  a  net  gain  in  precision  Cover  the 

A 

usual  estimator  R-y/x)  using  a  modified  estimator.  Such  as  estimator  is  based 


on  Double  sampling  (also  called  two  phase  sampling)  in  which  one  takes  the  ad¬ 


vantage  of  some  prior  knowledge  about  the  correlation  between  the  component 


-1 


characteristics  and  the  relative  costs  of  measurements  q  -  (q^  +  q^)  q^  and 


q  ■  qx(q^  +  q)  '  by  allowing  a  considerably  large  sample  on  the  cheaper  component 


(X  or  Y)  at  the  expense  of  a  small  reduction  in  the  size  of  the  joint  sample. 

The  idea  of  double  sampling  was  first  proposedby  Neyraan  (1938)  in  the  context  of 
stratified  sampling.  Since  then  many  authors  have  used  it  for  different  problems, 
see  for  example,  Robson  and  King  (1953),  Sukhatme  (1962),  Raj  (1963),  Sendransk 


[(1965,1967)],  Khan  and  Tripathi  (1967),  Sen  (1972). 

We  have  proposed  two  different  modified  estimators  of  the  ratio  R*Y/X  de¬ 
fending  of  whether  X  or  Y  is  cheaper  and  compared  them  with  the  usual  estimator 
* 

R*y/x  in  terms  of  bias  and  MSE.  Accordingly  two  different  cases  are  considered 
in  sections  2  and  3 


case  (i) :  q  -  q^  (q^  +  q^,)”1  >1  (X  cheaper) 

case  (li):  q*  »  qx(qQ  +  qy)”1  >  1  (Y  cheaper). 


In  section  3  the  method  of  approximation  of  Bias  and  MSE  has  been  discussed.  Thl 
is  in  the  same  spirit  as  in  David  and  Sukhatme  (1974)  except  that  asymptotic 
in  the  context  of  douhle  sampling  has  been  clarified  and  unified  with  the  usual 

A 

asymptotic  for  R. 


2.  THE  CASE:  q  -  qy(qQ +  qx)"’1  >  1 

As  mentioned  earlier,  the  modified  estimator,  proposed  in  this  section, 

A 

will  be  compared  with  the  usual  estimator  R»y/x  based  on  a  s.r.s.  (w.o.r.)  whose 


size  n  is  pre-determined  by  the  budgetary  constraint:  n-B(q_+q  +q  )_1. 


Our  modified  estimator  is  based  on  a  double  sampling  scheme  (or  a  two  phase 
scheme)  in  which  s.r.  sampling  (w.o.r.)  is  used  in  both  phases.  If  we  allow  a 
reduction  of  order  an,  ae  (0,1),  in  the  size  of  the  joint  sample  (that  we  use  for 

a 

R)  i.e.  if  we  allow  a  sacrifice  of  an  observations  on  the  costlier  component  Y, 
we  can  take  qan  additional  observations  on  X;  we  implicitly  assume  an  and  qan  to  be 
positive  integers.  We  shall  denote  by  n  and  n.  respectively  the  size  of  the 
smaller  (than  nl  sample  on  Y  and  the  size  of  larger  sample  on  X;  we  shall  frequently 
use  the  identities 


ng  ■  n-an  ■  n(l-a) 

n  »n+qan  ■  n(l+qa).  (2.1) 


Let  x.,x  ,...,x  be  the  sample  of  X  -  measurements  for  the  units  selected  in 

X  4  K 

the  first  phase  leading  to 

n„  n„ 


-  v2 


k.  ■  —  7  x . ,  and  s.  °  - r  7  (x  -xa) 

1  nl  j-1  ^  1 x  "i"1  j-1  3  1 


(2.2) 


which  are  known  to  be  unbiased  estimators  of  X  and  Sx  (Cochran  [1977]  pp.  19-27), 
where 


,  H  ,  ,  N 

5-*M  ““  sx ■  l <V»  • 

i-1  1  i»l 


(2.3) 


It  is  also  well  known  that 


E(v  *>2  -  u-f*>  4/n* 


(2.4) 


where  f^  ■  n^N  \ 

Let  {j^,j2»»*«»j  )  be  the  sample  of  units  drawn  from  the  set  of  first  phase 

s 

units,  j,  e  (1,2,. . .  ,n},  k«l,2,...,n  ,  and  let  their  Y  -  values  be  y.  ,y.  ,...,y.  , 
k  s  h  4  i3 


leading  to 


we 


Instead  of  the  Bias  B(R^)  •  E(R^-R)  and  the  mean  squared  error  MSE  (R^)  *  E  (R^-R) 

A 

will  consider  the  relative  bias  and  the  relative  mean  squared  error  of  R^  as 
well  as  of  any  other  estimator.  Denoting  them  respectively  by  RB  and  KMSE, 
we  have 

A 

*  R.  i 

RBCR^  -  EC-^-1)  "  j  BCR^ 

A 

RMSE(R  )  -  EC-~-l)2  -“MSECR.)  (2.9) 

1  R  IT  1 

St  A  A 

and  denoting  their  1  order  approximations  respectively  by  RB^ (R^)  and  RMSE^CR^) 

they  turn  out  to  be 

“i<v  • 

■  tiS+qa  "  (CX  ~  “ScSt* 
s  l 

“n  "  l+qa)CY  +  ^1+qa  "  f)CXY} 

where 

f  -  n|N 

cx  -  V*'  Sr  *  V5 

4  -  4cx-2ocxV  <2-io> 

The  errors  in  these  approximations  satisfy 

IRBCR^  -RB1(R1)|  -0(n"2) 


C2.ll) 


6 


IRMSEO^)  -BMSE1(R1)|  -  0(n"2). 


(2.12) 


A  proof  In  a  slightly  general  setting  and  with  some  mild  assumptions  appears  In 
section  3.  Asymptotics  for  the  context  of  double  sampling  has  been  clarified  and 

A 

unified  with  the  visual  asynptotics  for  R. 

A 

Comparison  with  R 

It  is  well  known  (see  Sukhatme  and  David,  1974)  that  a  first  order  approximation 

A  A  A 

of  the  relative  Bias  B(R)  and  the  relative  MSE(R)  of  R  are 


“i«>  '  £-5>(cx-pCxV 

-  ^Cl-f)(cJ-pCxCY)  (2.13) 

and 

““i®  •  i  - 11 4 

-i(l-f)c4y  .  (2.14) 


The  errors  satisfy 


and 


|RB(R)  -  BX(R)  |  -  0(~) 


RMSE  (R)  -  MSE1(R)  -  0  (rj)  . 

n 


(2.15) 


(2.16) 


Expressing  the  RMSE^R^)  in  terms  of  RMSE^R),  using  (2.10)  and  (2.14),  we 
obtain 


RMSE^)  -  RMSE1(R) 

CXY  ,  _c[a 


-  f  gCrl\  \ 

n  v  (1-a)  (1+qa)  a  q+qC^ 


{- 


(2.17) 


-  i  •— : 


Sf.' 

fell 


v.  •  .'•Vi  . ->  A.'< 


7 


where 


with 


CXY_CY  2 

C  2  "  ScIy  “  2pCx|y 

CX|Y  "  CX/CY  CY|x  "  CY/CX’ 


In  case  C  is  positive  and 


would  perform  better  than  R  as  long  as  a,  0  <a  < 1,  satisfies 


(2.18) 


(2.19) 


In  terms  of  the  correlation  coefficient  p,  Cx|y  *  Cx^CY  and  Cy|x  *  1 / cx | y 
conditions  (2.18)  and  C2.19)  can  be  expressed  as 

2pCy.«Cx-c£/qCx  (2.20) 

and 

a  < 

Observe  that  if  coefficient  of  variation  of  the  variable  X  is  twice  as  much  as  the 
coefficient  of  variation  of  the  variable  Y,  and  the  cost  ratio  q  is  moderately 
large  (say  q^.5),  the  condition  (2.20)  holds.  In  terms  of  p,  whenever  q  is  large 
and  p  <^0,  (2.20)  is  clearly  satisfied. 

Assuming  that  we  are  in  a  practical  situation  in  which  (2.18)  is  satisfied, 
the  best  choice  of  a  can  be  found  by  minimizing  (2.17)  w.r.t.  a.  Taking 
the  first  derivative  of  the  bracketted  expression  in  (2.17)  and  equating  it  to 


zero,  we  are  led  to  solve  the  quadratic  equation: 


(q-C)a2  +  2Cl+C)a  +  (l-qC)q_1  -  0. 


(2.21) 


The  roots  of  (2.21)  turn  out  to  be 

-2(1+C)±  A 
2  (q~C) 

where 


(2.22) 


A  ■  2j^C  (1+q)/*^  . 


Denoting  the  critical  points  By  Oq  and  a^,  we  get 


«  SaSc-  1  _  Sc-  (Sq)"1 

q+v^qy^f  Sc+Sq 


n  .  S~ a^f+1 
1 


(2.23) 


Due  to  the  assumption  q>l  and  qC  >  1  (vide  (2.181),  it  follows  that  a^  i  (0,1) 
otg  e  (0,1).  The  second  derivative  of  (2.17)  with  respect  to  a,  evaluated  at  a^, 
turns  out  to  be 


— ~(JKMSE-  (R.)) 
3a  A  A 


2q2<$ 


n(l+qaQ)' 


{C  +  t 


1+qV3  _1_, 

l-a_  ;  2* 

0  q 


(2.24) 


which  is  positive  for  C>_0. 

Hence,  under  the  condition  (2.18),  the  optimum  size  of  the  joint  sample  is 


n_  “  nU-aJ 


(2.25) 


and  the  size  of  the  additional  sample  on  X  is 


nq  Q  =  n 


Sq  Jc-i 


l  +  /C/q 


(2.26) 


The  RMSE^(R^)  corresponding  to  the  optimum  choice  of  a  turns  out  to  be 


[RMSE  (R  )  ]  =  RMSE  (R)  -  - — -  (/C )2. 

1  U  1  _  /I  .  J-n  JI 


(2.27) 


n(l  +  — ) 

q 


-1 


Ignoring  the  terms  of  order  N  , 


1x2 


W“l(Y  ]0  *  R  n  n 


2  CS  _  S 


(1+-) 

q 


(2.28) 


and  when  q  is  very  large 


(2.29) 


nX 


implying  that  for  a  small  reduction  in  the  joint  sample  (see  (2.25))  a  very 
large  sample  of  additional  observations  (see  (2.26))  enables  the  modified  esti¬ 


mator  behave  as  if  its  denominator  were  the  population  mean  X.  In  terms 


_ 2  2 

of  the  parameters  R  =  Y/X,  S^,  and  p,  the  optimum  [MSE^(R^)]q  is 


‘“Wo  '  n(l!,)x"  Sy  Vr2S*-2*>S_S„  V 


X  Y 


-t4 


(2.30) 


The  relative  efficiency  of  the  optimum  modified  estimator  R^  (say  (R^)q  w.r.t.  R 
is 


10 


R.E.  CRj^q 


MSE^R) 


rmse^Cr) 


[mse^r^  [rmse1(r1)]0 


(l-Hq)[l  +  Cx|Y-2  Cx|y^  "n*1*  CxlY~2 
[^■+/CX|Y  "  2pCx|Y  '  “  N  ^  +  Cx|Y~  2pCx|Y^1+q) 


(2.31) 


In  situations  in  which  *  0,  we  have 

N 


[MSE^R^  -  MSE1(R1)0 


1  r2 


n(l+q) 


X  [v£  SY  yRZS^-2R  SxSy  ] 


(2.32) 


“‘Vo  * 


C1-N)[1-^Cx|y-2  Cx[y] 
+\/cx1y“  2pCx|y  1 


(2.33) 


Further  if  cx  Y  *  Cx  Cy  ■  1 


:CR1)  * 


]2 


(2.34) 


3.  APPROXIMATION  OF  BIAS  AND  MSE 

a  A  ^  mm  " 

Let  X  and  Y  denote,  respectively,  unbiased  estimators  of  X  and  Y,  based  on 
a  finite  number  of  randomly  chosen  units.  We  assume  taht 


for  all  possible  realizations  of  the  sampling  scheme,  and  we  propose  the  ratio 


R  -  Y/X 


(3.1) 


as  the  estimator  of  the  population  ratio 


V 

O 


A'- 

A-- 


m 


1 


r.  * 
Co 


m 

P 

I 


k 


m 


V  w 


R  -  Y/X. 


The  bias  of  this  estimator,  in  terms  of  the  random  relative  deviations 


5  -  (Y  -  Y)/Y 


and 


€  -  (x  -  x)/x 


is 


-1, 


BCR)  “  E(R-R)  «RE[(6-e)(l+e)  ]. 


th 


If  we  define  a  k  order  approximation  of  B(R)  as 


2k 


Bk(R)  -  RE [ (6-g )  l  (-e)1"1], 
i-1 


then  the  error  of  this  approximation  turns  out  to  be 


efcB  -  B(R)  -  BfcCR) 


2k, 


Y  E[(l/X)C6-€)€  ]. 


The  mean  squared  error  of  R,MSE(R),  in  terms  of  6  and  e  is 

,2 


MSE(R)  -  E(R-R) 


R2E[(6-6)2(l+€)"2]. 


th 


If  we  define  a  k  order  approximation  of  MSE(R)  as 

2R-1 


MSE  CR)  -  R^E [ (6-e ) L  \  iC-e)1-1] 

K  i-1 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


C3.6) 


(3.7) 


(3«81 
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then  the  error  of  this  approximation  turns  out  to  be 


ekMSE  -  MSE(R)  -MSE^R) 


-Y  E[a/X)2(5-6)2C2k€2k“1+(2k-l)€2k)]  .  (3.9) 


In  the  sequel  we  give  bounds  for  e^B  and  e^MSE  for  each  of  the  estimators 
considered  in  this  report.  In  doing  so  we  follow  the  techniques  and  arguments 
of  David  and  Sukhatme  [1974].  We  also  make  use  of  the  following  order  relation 
the  proof  of  which  can  be  found  in  the  work  of  these  authors;  for  non-negative 


integers  r  and  sf 


Y1^3  E 


«V>  -  f 


Otl/n^^2}  if  t+s  is  even 


0{l/n(t+8+1)/2}  if  rfs  is  odd 


(3.10) 


where  the  random  'relative'  deviations 


cx-xyx  -  (x-x)/x , 


(Y-Y)/Y  -  (y-Y)/Y, 


are  based  on  the  sample  means  x  and  y  of  a  simple  random  sample  (without  replacement) 

*  A 

of  size  n.  However,  when  the  sample  means  X  and  Y  are  based  on  a  two-phase  sampling 
design  in  which  a  large  sample  of  n^  units  are  drawn  in  the  first  phase  to  com¬ 
pute  one  of  the  means  and  a  subs ample  of  ng  units  are  drawn  out  of  first  phase 
units  to  compute  the  other,  a  similar  order  relation  follows. 


REMARK 


_1.  In  the  case  of  Rj^  *  y^x^. 


5  -  «.  -  G  -Y)/Y  -  E(6_ In. ) 


where 


5s  ■  cV*)/?  "  +uVV/f} 

and  EC* | n^)  denotes  conditional  expectation  given  the  Y  -  values  of  the  first  phase 
units. 

REMARK  2.  In  the  case  of  R  *  y  /x  , 

2  x,  s 

€  -  Cx  -Xi/X  -  EC€s|n£)  (3.12) 

where 

6S  "  CV5)/S  "  +{(vV/5}- 

*  A 

We  now  study  the  order  relation  for  the  estimator  R^  The  same  for  R2  can  be 
studied  similarly.  For  the  positive  and  even  integers  r  and  s  the  computation  of 
the  order  of  E(yg-Y)r(xjl-X)S  proceeds  as  follows; 

CE[6^®]  -  E[(y8-Y)2a(xrX)2b] 

with  r  *  2a4s  -  2b,  where  a  and  b  are  non-negative  integers  and  C  -  XsYr.  Let 


and 


We  note  that 


X»_x*  Ay. 


y„  -  ye 


aV 

s  l 


Wx1)2b(iJij|i+ iy£)2a, 


(3.13). 


which  in  turn  can  be  expressed  as 


14 


rs  /A—  N2br  r  *2aw  —  \2I/A-  x2a-2i 

,u  -  C4V  i  I  (iyt> 


.  %  ,2a  .  -  ,21+1,,-  ,  2a-2i-l, 

+  1i0<21+l)W>'s|l)  <4V  1 


(3.14) 


Now  taking  the  expectation,  conditionally  given  the  first  phase  units  and  adopting 
the  order  relation  (3.10)  also  for  this  situation,  dispensing  momentarily  with  the 
justification  for  doing  so,  we  have 


,  -  .2b r  “  /2a.  1  .  21,  -  ,2a-2i 

(AX  )  [  l  C  )  0(— )  (Ay  A 

*  l-o  Z1  AT  * 


1  .21+2,.-  .2a- (21+1) 


+  IA1+ 1» 

1*0  /  n 


(3.15) 


Next,  the  expectation  with  respect  to  the  first  phase  units  yields  the  unconditional 
expectation  as 


f  (£>  o(-M21  o^)2at2b-21 

L-0  JiT  ST 

a  l 


+  ?<£*>  0(-i=)21+2  otp 

1"0  /n  /n. 

s  1 


1  .21+2  1  . 2a+2b-(2i-l)+l 


which  in  turn  is  of  order  0(l/«^T)^a+^> ,  yielding 


C  E[  Sr€3]  ■  0(— )I+S,  r+s  even; 

3  1  Sn 

since  n  •  (l-a)n  and  n.  -  (l-qa)n  with  0<  a<  1.  Similarly  we  get  the  order  relation 

S  X, 

for  the  other  cases  of  r+s.  The  justification  for  these  order  relations  are  as 
follows . 

A  A 

In  order  to  compare  R  with  R^,  their  respective  settings  for  asymptotics  must 
be  interlinked  and  this  In  both  situations  1-  1,2.  For  doing  so  we  implicitly  follow 

A 

the  setting  of  David  and  Sukhatme  [1974]  for  the  asymptotic  study  of  R  according  to 


which  a  sequence  of  finite  populations  S.  ,S0,..,,S  of  sizes  N, ,N„,...,N  is  drawn 

i  i  n  x  l  n 


from  a  bivariate  super  population  S  -  { (u^ »v1) »^u2’v2^ » • • • > •  The  sequence  of  sizes 

is  strictly  increasing  so  that  n-  t^N^,  lim  t  «  t,  0<t<l  and  0<  t  <1  for  every  n 

n-*- 00 

The  sampling  scheme  adopted  in  the  case  of  R^  (i «  1,2)  is  two  phase  sampling. 

In  the  first  phase,  a  simple  random  sample  (w.o.r.)  of  size  n  -  (l+qa)n,  qa>  1, 

is  drawn  from  S  .  We  assume  that  n„ <  N  for  every  n,  so  that 
n  l—  n  } 

n„  »  (l+qo)t  N  ■  t  N 
l  ^  n  n  n  n 

*  *  * 

The  condition  n-*-  *  implies  n^-*-®  and  also  that  t  -►  t  *  (l+qa)t  with  0<  t  <1  and 
* 

0  <  t  <1  for  every  n.  In  the  second  phase,  a  simple  random  sample  (without 
replacement)  of  size 

l“Qt 

n  -  (l-a)n  ■  — —  n  ■  X  n. ,  0  <  a  <  1, 

s  1+qs  ini*  * 

X  »  X  ■  (l-a)/l+qs,  is  drawn  from  a  sub-population  consisting  of  n  units, 
n  i 

The  sizes  of  these  sub-populations  are  strictly  increasing.  The  condition  n  -►  ® 

s 

is  implied  by  the  condition  n-*- «  and  we  also  have  lim  X  *  X  with  0<  X-  X  <1 

n  n  1 

for  every  n.  The  settings  for  asymptotics  being  linked  in  this  manner  the  content 
of  Sukhatme  and  David  [1974]  applies  in  verbatim. 

Bias  and  MSE  of 

Following  (3.4),  (3.6),  and  (3.7)  and  the  Remark  1, 

B(RX)  -REK^-e^O^)"1],  (3.16) 

i-1 

Bk(R1)  -RE[(6a-e4)  ]*  (3*17) 

and 

e^B  -  Y  EKl/^Ci^^f.  (3.18) 


' 

i 


As  mentioned  earlier,  we  assume  that  X-measurements  are  positive  for  all  the  units 


of  the  population.  Let  Xq  be  the  lower  bound  for  x^,  we  have  then 

e^I-  [E(6g€f)  -  E(e/k+1]  (3.19) 

and  using  the  order  relation, 

jekB|  <  0(l/n)k.  (3.20) 

This  bound  over  the  error  of  paproximation  of  order  k  shows  that  it  suffices  to 

A 

approximate  B(R^)  by  B^(R^)  for  sufficiently  large  sample  size;  letting  k«  1  in 
(  3.5  >.  we  get 

B^)  -  RE  [<V€Z)( l"€i)] 

■“iV^vVi1  (3-21) 

where 


Similarly 


and 


EC6s)  -  EEC6g|njl)  -  E(6 £) 


•  r  E(y  -?)-  0. 


E(S  £,)  -  ±  E[(y  -Y)(x  -X)] 


-  p  C  C 
x  y 


E(«2)  -  Jj  E(J  -S)2 

X 


-£-K 

i 


I1 

\\ 

4* 


we  have  therefore 


V\l8)  -  ^-?KCx-“CxV' 


As  to  the 


MSEC^) 


we  have 


MSECRX)  -  R2E[C«s-€t)2(l+e4)"2] 


MSE^)  =  R2E[<$g-ejt)2  l  i(-€)i-1] 


efcMSE  -  -  Y2  E[k/52)(5  -£  )2C2ke2k_1  +  (2k-l)e2k)J 


Noting  that 


l«kMSE|<-Sj|2k  E{<^-2«  «  -H?).**1) 


+  (2k- 1)  «Uq>2-ivX)*?,l 


we  observe  that 


lejMSEl  <0(l/n)] 


* 

and  consequently  approximate  MSE(R^ | R)  by  MSE^ (R^ | R)  which  turns  out  to  be 

MSE^RjjR)  -  E[C«s-€a)2] 


-  -  b C2  +  (—  -  b(C2-  2pC  C  ) . 


Computations  similar  to  those  involved  in  the  case  of  lead  to  the  following 


|ekB|  -  0(l/n)  , 


|ekMSE|  -  0(l/n)  , 


s  l  1 


MSEi(R2|R)  -  (f  -  ±>c*  +  (j-  -  ±)<cj-  2p  yy. 

8  At 


(3.29) 


(3.30) 


(3.31) 


(3.32) 
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